Conformal Invariance in Inverse Turbulent Cascades 



O 
(N 

o 

Q 



G. Falkovich^ and S. Musacchio^ 
^ Weizmann Institute of Science, Rehovot, Israel 
^ SCNRS, Lab. J. A. Dieudonne UMR 6621, 
Pare Valrose, 06108 Nice (France) 
(Dated: December 20, 2010) 

We study statistical properties of turbulent inverse cascades in a class of nonlinear models describ- 
ing a scalar field transported by a two-dimensional incompressible flow. The class is characterized by 
a linear relation between the transported field and the velocity, and include several cases of physical 
interest, such as Navier-Stokes, surface quasi-geostrophic and Charney-Hasegawa-Mima equations. 
We find that some statistical properties of the inverse turbulent cascades in such systems are con- 
formal invariant. In particular, the zero-isolines of the scalar field are statistically equivalent to 
conformal invariant curves within the resolution of our numerics. We show that the choice of the 
conformal class is determined by the properties of a transporting velocity rather than those of a 
transported field and discover a phase transition when the velocity turns from a large-scale field to 
a small-scale one. 
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Exceptional role of conformal invariance in theoreti- 
cal physics stems from the fact that most of non-trivial 
exact solutions of dynamical and statistical models can 
be traced to the existence of this symmetry. One of the 
most remarkable recent advances in mathematics was the 
discovery of Schramm-Loewner Evolution (SLE) and of 
the bridges it builds between different branches of physics 
[ij-Q. SLE is a class of fractal random curves that can 
be mapped into a one-dimensional Brownian walk and 
thus have conformal invariant statistics. SLE curves ap- 
pear at two-dimensional (2d) critical phenomena as clus- 
ter boundaries, thus revealing a statistical geometry of 
conformal field theories. 

In equilibrium, the statistical weight of a state does 
not depend on how the state was created. If a system is 
driven away from equilibrium by an external force then 
the probability of a given configuration depends gener- 
ally on the history and on the statistics of the driv- 
ing force. Turbulence statistics are thus generally force- 
dependent. All the more surprising was then the ex- 
perimental discovery that the isolines of vorticity in the 
2d Navier-Stokes turbulence and of temperature in the 
Surface Quasi-Geostrophic (SQG) turbulence belong to 
SLE 0, Q. That means that at least a part of turbu- 
lence statistics could be described in terms of a confor- 
mal field theory like equilibrium critical phenomena. In 
particular, nodal vorticity lines happen to be equivalent 
to the boundaries of percolation clusters Q, while the 
iso-temperature lines in SQG are equivalent to the do- 
main walls of S0(2) model (that of a 2d gaussian free 
field) 0. Having only two examples leaves wide open 
possibilities for different interpretations and hypothesis, 
particularly trying to relate the scaling and properties 
of the bulk field to the choice of the curve class for its 
isolines Q. Here we study additional models from the 
family and show that the SLE class is actually sensitive 
to the type of dynamics (i.e. velocity) rather than to the 
type of a field that is carried; that sensitivity manifested 



dramatically by the phase transition we discover. 

The class of models we investigate has been introduced 
in 0, [I]- It describes the evolution of a scalar field 9 
transported by an incompressible two-dimensional veloc- 
ity u ~ {dy'tp, —dx'ip), expressed via the stream function 
tp. The scalar field 9 is "active" because it is linearly 
related to V' and u. In Fourier space the relation reads: 
6{k) = \k\^"'ip{k). The system is thus governed by the 
equation 



dtO + {uV)9 = dt9 + {9,i^}^F + D 



(1) 



where {9,^p} = O^ipy — 9y^x, F and D are external forc- 
ing and dissipation respectively. Different values of m 
give different well-known hydrodynamic equations. For 
m — 2 one obtains two-dimensional Navier-Stokes (NS) 
equation, 6 being the vorticity. For m = 1 the field 9 
represents the temperature in SQG turbulence. Finally, 
for m = — 2 the model corresponds to that derived by 
Charney and Oboukhov for waves in rotating fluids and 
by Hasegawa and Mima for drift waves in magnetized 
plasma in the limit of vanishing Rossby radius (ion Lar- 
mor radius for plasma physics). 

At all values of m equation ((T|) possesses two positive- 
definite invariants for F = D = 0, namely Z — J 9^dx 
and E — J 6'ipdx/2. When the system is forced by an 
external source of scalar fluctuations F, with a correla- 
tion length £f ^ ^/kf, the existence of two conserved 
quantities causes double turbulent cascade. The sign of 
TO determines the direction of the cascades. For m > 
the "energy" E is transferred toward large scales i > (f 
giving rise to an inverse cascade, and the "enstrophy" Z 
flows toward small scales. The cascades are reversed for 
TO < 0. 

Here we focus on the range of scales corresponding 
to an inverse cascade. Dimensional argument based on 
the assumption of scale-independence of the flux of en- 
ergy in the inverse cascade (for to > 0) gives the scal- 
ing exponent h = (2 — 2to)/3 for the the increments 
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FIG. 1. Scaling exponents h of the scalar field Q (circles) 
and stream function (triangles), and universality class k 
(squares) for various values of m. 



FIG. 2. Pdfs of scalar increments 5Q(r) at different distances 
r within the inertial range for m = 3/2, m = 1/2, m — > 
(from top to bottom) 



6rQ = 6i(a; + r) - d{x) - r'' and /i = (2 + m)/3 for 
the stream function V'- For m < similar argument 
gives the scaling exponent h = {2 — m)/3 for the field 
9, and h = {2 + 2m)/3 for the stream function ip. Our 
direct numerical simulations support these predictions. 
The scaling exponents of the fields 9 and ip are shown in 
Figure [Has a function of m. When m ^ —m, fields 9 and 
tp exchange their scaling exponents (and the respective 
cascades change directions). 

The first remarkable discovery of conformal invariance 
in turbulence has been made for the zero-vorticity lines 
in Navier-Stokes turbulence i.e. for m = 2 [4|. Zero- 
vorticity regions correspond to Aijj = 0. i.e. to a har- 
monic stream-function and are invariant with respect to 
conformal transformations (which thus map streamlines 
into themselves). One may think that conformal invari- 
ance of zero-vorticity lines is a remnant of the invariance 
of zero-vorticity domains and is peculiar for m = 2. How- 
ever, conformal invariance of the isolines was then dis- 
covered for m = 1 [3], where one does not recognize an 
analogous property of zero-9 domains. It is then tempt- 
ing to relate conformal invariance to the properties of 9 
which are common for all m. The main property seems 
to be the fact that 6* is a Lagrangian invariant of the flow 
and determines the symplectic structure. For example, 
if we denote R ~ {X, Y) the initial (Lagrangian) coor- 
dinates of the fluid particles then the extremum of the 
action I = J S{t) dt with 

S'2(r) = J 9{Ii)x{R,t)y{R,t)dR 

-i y"0(Ri)0(R2)hi|r(Ri)-r(R2)|dRidR2 (2) 

gives X = dyip and y — —dxip which is equivalent to ^ 



for m — 2. Generally, 

S„,ir) = J 9iR)xiTl,t)y{Tl,t) dR 

y"0(Ri)0(R2)|r(Ri)-r(R2)|"-'dRidR2 . (3) 

In other words, the energy E = j '4){y)9{y) dr /2 — 
J 6l(ri)6l(r2)|r(Ri) - r(R2)r-2 rfridr2/2 is the Hamil- 
tonian. It is tempting to conjecture [9j that zero-9 lines 
are special since the Hamiltonian description is singular 
(non-invertible) there. However, at negative m, 9 is a 
large-scale field and its isolines are not fractal (have di- 
mensionality 1). It is then natural to study the properties 
of the isolines of ip which are fractal now. We show be- 
low that at TO = —2 the isolines of ip seem to have the 
same statistical properties as the isolines of at m = 2, 
despite the fact that tp is not a Lagrangian invariant and 
the equation has no symmetry m ^ —m. 

To investigate the statistical properties of the scalar 
field 9 we solved numerically eq. [T] on a doubly peri- 
odic square domain of size L = 27r at different resolution 
N'^ = 1024^, 2048^. The scalar fluctuations are sustained 
by a Gaussian, (5-correlated in time, random forcing, 
peaked around wavenumber fc/ = 100. Dealiasing cut- 
off is set to kt = N/3. Time evolution was computed by 
means of a second-order Runge-Kutta scheme, with im- 
plicit handling of the linear dissipative terms. The direct 
cascade of enstrophy is halted at wavenumbers fc > fc/ by 
means of a hyper-viscous damping (— l)P~^i'p_i V^^ of 
order p = 8. Statistically steady state in the inverse cas- 
cade is obtained by removing the energy at large scales 
with a linear friction term —ri9. Note that for to > the 
characteristic times of the inverse cascade process scales 
as Ti ^ ^(4-m)/3 j -[-jj^g cascade slows down as m goes to 
zero. This phenomenon limits the resolution achievable 
in numerical simulations. 
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FIG. 3. Power spectrum of the scalar field 6 for m = 1/2 
Dashed line is the slope k~^^^ . The inset shows the energy 
flux. 
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FIG. 4. Compensated power spectrum of the scalar field 9 for 
m ^ 0. Here m — 10~^, 10~^, 10"'^. In the inset it is shown 
the power law behavior of the amplitude C(m) ~ rn~^/^. 



The scalar field resulting from numerical simulations 
with < TO < 2 is scale invariant, as confirmed by the 
perfect collapse of the probability distribution functions 
(pdfs) of scalar increments (5^6* for different r, see e.g. 
Figure [H Note that the pdfs are non-Gaussian. The 
scalar field also has a power-law spectrum for k < kf in 
agreement with the prediction: 



(4) 



where C = (4to — 7)/3, and e is the flux of energy (see 
e.g. Figure [3]). 

The limit to — > of the active scalar model is singular. 
Indeed for to — the two fields 9 and tp coincide, and 
the advection term {9, ip} in eq.[l] vanishes. Therefore no 
turbulent state can be produced and the field 9 is sim- 
ply determined by local balance between forcing and the 
dissipation at exactly to = 0. Conversely, for arbitrary 
small values of m we find a turbulent cascade with power 
law spectrum in agreement with eq. |3](see Figure S]). As 
the parameter m goes to zero, the amplitude of the scalar 
field diverges, to compensate for the less efficient transfer 
of energy in the cascade. This is signalled by the power 
law behavior of the analogous of Kolmogorov's constant 
for the spectrum C{m) ^ m~'^/^ (see inset of Figure 

To study the limit m -> let us write the l.h.s. of ([1]) 
in fc-space, and use the symmetry o fc — j: 



d9k 
mdt 



= ^T.^k,j]r^9,9u-i 



(5) 



limit TO — > 0, equation ([5]) has still the form of a transport 
equation with the link between 9 and the stream function 
being ip{k) = — ln|fe|6'(fe). Renormalizing t rat one 
gets 



d9k 
dt 



\ Y}k,j] H3l\k - j\)9j9k-j +F + D 



(6) 



Numerical integration of eq. ^ produces an inverse tur- 
bulent cascade with power law spectrum Pe{k) ~ k~'^ 1'^ . 
(see Figure [5l) In the range of scales of the inverse 
cascade the field 9 is self similar with scaling exponent 
h — 2/3, as confirmed by the re-scaling of the pdfs of 
scalar increments (see Figure [21). 




where [fc,j] = fcij2 — ^2^1, = |fc| and j — In the 



FIG. 5. Power spectrum of the scalar field Q for the asymp- 
totic model rn ^ 0. Dashed line is the slope k~''^^ . The inset 
shows the energy flux. 

Numerical investigation of Navier-Stokes (NS) equa- 
tion 4] and Surface Quasi Geostrophyc (SQG) model Q 
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FIG. 6. Perimeter P of zero-isolines versus gyration radius 
L for m = 1/2. Forcing length scale is r/ ~ 0.06. In the 
inset we show the local slope of perimeter P before and after 
randomization of the phases of the scalar field 9 (solid and 
dashed line respectively). 

have shown that for two pecuHar cases, namely for m — 
2, 1 the zero-isohnes of the scalar field are statistically 
equivalent to SLE i.e. could be mapped to Id Brownian 
walk. The class SLE is characterized by the respective 
dimensionless diffusivity k [H, 0]- In particular for NS 
the zero-vorticity isolines belong to the same universal- 
ity class of critical percolation and are equivalent to SLE 
curves with k = 6. For SQG the zero-temperature iso- 
lines are SLE curves with k = 4. It is therefore natural 
to ask if the properties of conformal invariance for the 
zero-isolines is a general property that holds for arbitrary 
values of m. 

To investigate this issue we consider the connected re- 
gions of positive/negative sign of 6, The boundaries of 
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FIG. 7. Perimeter P of zero-isolines versus gyration radius L 
for m ^ 0. Forcing length scale is r/ ~ 0.06. 




FIG. 8. Statistics of the driving ^(t) for m = 1/2 



these clusters are closed loops formed by the zeio-0 iso- 
lines. 

For < m < 1 the scalar is a self-similar rough 
field with scaling exponent < ft. < 1. The relation 
between the scaling exponent h of a height function and 
the fractal dimension D of its isolines was suggested in 
[13]: (3- ft)/2 = {7+2m)/6 for < /i< 1 and < m < 
1. Let us stress that this is not the fractal dimension 
of the iso-set (known to be equal to 2 — /i for /i > 
and to 2 for /i < 0) but that of a single long isoline. 
One can thus conjecture the relation k = 4(1 + 2m)/3. 
Indeed for m = 1 it was found Q that the zero-isolines 
are SLE curves with k = 4, in agreement with the above 
prediction. Nevertheless, in Figures [6] and [7] we show 
the fractal dimension of the isolines for to = 1/2 and for 
the asymptotic model to — > 0. It both cases the fractal 
dimension measured is not agreement with the prediction 
D = {7 + 2m)/6 but is compatible with D = 3/2. 




FIG. 9. Statistics of the driving ^(t) for m ^ 
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FIG. 10. Variance of winding angle as a function of the 
lenght I along the isoline for m = 1/2 



FIG. 12. Variance of winding angle as a function of the 
lenght / along the isohne for m = 1/2 (squares), m — 3/2 
(circles), m = 2 (triangles). 



Following the procedure described in j5| from the zero- 
field lines we obtain an ensemble of curves in the the half 
plane which are expected to converge in the scaling limit 
to chordal SLE. Then we extract the driving ^(t) of the 
corresponding Loewner equation. As shown in Figures [5] 
and [3] the driving has Gaussian statistics with variance 
(f^(i)) ^ Kt. For all the cases considered with < m < 1 
we found k — 4, which is in agreement with the fractal 
dimension observed. 

As a further test we study the statistics of the wind- 
ing angle of the zero-isolines. The winding angle (p is 
defined as the degree with which the curve wind in the 
complex plane about a point w <p{z) — arg(z — w)^l, 12[. 
The asymptotic distribution for the winding angle at long 
distance £ along the curve is Gaussian, with variance 
k/(4 -I- k/2) log(^). As shown in Figures [10] and [TT] we 
found that its variance grows like 2/31og(^), thus sup- 
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FIG. 11. Variance of winding angle as a function of the 
lenght I along the isoline for m — > 



porting the conjecture k = 4, and is not compatible with 
the prediction k = (4 -I- 8m)/3. 

This model provides an example of non-trivial relation 
between the scaling exponent of the field 9 and the frac- 
tal dimension of its isolines. For < m < 1 the scaling 
exponent varies in the range < h < 2/3, but the fractal 
dimension remains constant D = 3/2, at variance with 
what one would expect from the relation D = {?> — h)/2 
which holds for Gaussian random field. The crucial dif- 
ference is that the scalar field 9 is not random, but is the 
result of turbulent dynamics. 

Our findings can be understood in Lagrangian terms. 
The scaling exponent of the velocity field is hy = (m — 
l)/3. For m > 1 the scaling exponent is positive, and 
therefore velocity difference scales as Sv{i) oc 
Two Lagrangian trajectories moving in such velocity field 
separate according to the Richardson law i{t) ^ t^^ (4-m) ^ 
Conversely for < m < 1, the exponent is negative (that 
is the velocity is a small-scale field like vorticity in Navier- 
Stokes), and velocity differences are independent of the 
separation Sv{£) ~ Vrms- Lagrangian trajectories will 
separate as £(t) (x t. Perimeter P and gyration radius 
L of clusters can be related by assuming that their ratio 
P/L, which is proportional to the number of folds, grows 
as a random walk, i.e. as t^^^. Gyration radius grows as 
two-point distance L{t) oc £{t), which gives P cx Lt^^'^ cx 
^(io-m)/6 foj. m > 1 and P oc £3/2 for < m < 1. 

The property of conformal invariance of the isolines is 
therefore determined by the underlying dynamics of the 
field. As a test we took the field 9 and randomize its 
phases in Fourier space. This procedure does not change 
the scaling exponent of the field, but destroys all the 
correlations generated by the turbulent dynamics. The 
isolines of this randomized field are no more conformal 
invariant, but their fractal dimension recover the predic- 



6 




FIG. 13. Statistics of the driving ^{t) for the stream-function 
isohnes for m = — 2 

tion for Gaussian random field (see inset of Figure [6]). 

Under the hypothesis that the isoHne of the scalar field 
9 are SLE curves one can obtain a conjecture for their 
universality class k from the formula for the fractal di- 
mension = 1 + 2/k oi the outer perimeter P, which 
holds for K > 4 One obtains k = 4 for < m < 1 
and K = 12/(4 — m) for m > 1, which are in agreement 
with our findings (m = 3/2,1/2 and m — >■ 0) and with 



previous results (m = 1,2). 

Approach based on Schramm-Loewner Evolution pro- 
vides a refreshingly novel geometric insight into the 
statistics of turbulence and hints at deep symmetry as- 
pects of 2d flows which we are yet far from understanding. 
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